Interaction of nonuniform elastic waves with two-dimensional electrons in 

AlGaAs — GaAs — AlGaAs heterostructures 
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An interaction of a double layer electron system realized in an AlGaAs — GaAs — AlGaAs 
heterostructure with nonuniform elastic modes localized in the GaAs layer is considered. The 
dependence of the coupling constant on the ratio between the thickness of the GaAs layer and the 
wavelength is calculated for the wave vector directed along the [110] axis, the displacement vector 
lying in the (lIO) plane and the interface boundaries parallel to the (001) plane. It is shown, that 
the coupling constant reaches the maximal value at the wavelength which is of order of the thickness 
of the GaAs layer. The renormalization of the velocity of the elastic modes is found for the case, 
when the electron system is in the fractional quantum Hall regime. It is shown, that for certain 
modes the dependence of the velocity shift on the wave vector is modified qualitatively under a 
transition of the electron system to the state, which corresponds to the Halperin wave function. 

PACS numbers: 71.10 Pm, 73.40.Hm, 74.20. Dx 



Recently the surface acoustic waves (SAW) have been 
widely used for an investigation of the dynamical prop- 
erties of two-dimensional electron layers in heterostruc- 
tures AlGaAs |l|,g]. Since the AlGaAs is a piezoelectric 
material, the SAW generates an alternating electric field. 
The interaction of that field with the two-dimensional 
electrons results in a renormalization of the velocity and 
the absorption of the SAW. The frequency and the mo- 
mentum dependence of the conductivity of the electron 
system can be obtained from the measurements of the ve- 
locity shift and the absorption coefficient. In particular, 
the power of the method has been demonstrated within 
the study of the quantum Hall systems for which the con- 
ductivity depends considerably on the external magnetic 
field. 

Since the electron layer is placed at a certain dis- 
tance do from the surface of the sample, the coupling 
of the SAW with electrons depends on the parameter qd^ 
(where q is the wave vector of the elastic mode) . An expo- 
nential reduction of the matrix elements of the coupling 
takes place if the parameter qd^ increases. Therefore, the 
SAW method is not applicable for a study of the dynam- 
ical properties of the electron system at large q. In that 
case a nonuniform elastic wave, which is localized near 
the electron layer due to the acoustic nonhomogeneity of 
the heterostructure, can be used. 

To specify the system, in which the elastic wave may be 
localized, we consider the heterostructure AlxGai^^As — 
GaAs— AlxGai^xAs which incorporates two electron lay- 
ers situated at the AlGaAs — GaAs interfaces. Similar 
heterostructures have been used for the investigation of 
the quantum Hall effect in the double-layer electron sys- 
terns li. 

Several types of the elastic modes localized in a central 
layer exist in the structure considered. The interaction of 
the two-dimensional electrons with the nonuniform trans- 
verse elastic wave with the displacement vector parallel 



to the interface boundaries has been studied in p|. In 
this paper we continue the investigation of Ref. [^ and 
consider the nonuniform elastic waves for which the dis- 
placement vector is in the plane determined by the wave 
vector and the normal direction to the interface bound- 
aries (these waves as well as the SAW are polarized el- 
liptically) . We found that the interaction of such waves 
with electrons may be much stronger then of the waves 
considered in Ref. B. 



I. GEOMETRY OF THE MODEL AND THE 
DISPERSION EQUATIONS 

Let us consider the system for which the GaAs 
layer with the thickness 2a is situated between two 
AZo.sGao.r^s layers with the thicknesses much larger 
than the wavelength of the elastic mode. We specify the 
case of the interface boundaries parallel to the (001) plane 
and the wave vector of the elastic mode directed along the 
[110] axis. We will use the fact that the elastic moduli are 
practically independent on the concentration of Al (we 
set them the same for the both media). The acoustical 
nonuniformity of the system is caused by the difference 
of the densities (/9i=5.3 g/cm'^ for GaAs, p2=4.8 g/cm'^ 
for A^o.sGao.T^s). 

We consider the elastic mode with the displacement 
vector components 



Ui{v,z,t) = Ui{z)e 



ic\v — iUqt 



(1) 



the 



where ojq = vq, v is the elastic mode velocity, i = x, z, 
X axis is chosen along the [110] direction, the z axis, along 
the [001] direction, the r vector is in the (001) plane. The 
wave equations for the Ui(z) components have the form 

(C449^ - C^iq^ + PaUJ^)Uj; + iq{ci2 + Cn)dzUz = 



iq{ci2 + C44:)dzU.j; + [ciidl 



caq 



p^u?)u,_ = 0, (2) 



where Cn — 0.5(cii + ci2 + 2C44), a— 1,2 corresponds to 
the medium number. The sohition of Eqs.(0) has the 
form 



'(^) = Xl^ScexpyfeZg, 



where y^ are the roots of the equation 



y'^ + 26„y2 _^ ^^ 







(3) 



(4) 



with 



ba = 



1 



2cilC44 



[(Cl2 + 044)^ + ClliPaV^ ~ Cii) 
+ C44{PaV'^ - C44)] 
-{PaV^ - Cii)iPaV^ - C44). 



C11C44 



(5) 



If we take into account the equivalence of the elastic 
moduli for two media, the boundary conditions are re- 
duced to the requirement of the continuity of Ui and dzUi 
at the interfaces The localized mode corresponds to the 
solution for which the displacement approaches to zero 
at z ^ ±c». The structure of the localized solution in 
the AlGaAs layers is similar to the structure of the SAW 
on the surface of the cubic crystal (see, for instance, Ref. 
|g[ ) . Such a solution arrears if Eq. (|j) at a = 2 does not 
have imaginary roots. If the elastic moduli satisfy the 
inequality 



(Ci2 + C44)^ - Cii(Cii - C44) < 0, 



(6) 



the localized solution emerges at w < c^a/ p2- If the op- 
posite inequality is satisfied, the velocity of the localized 
mode V < Vm2, where Vm2 is the root of the equation 
D2(v) —h\~C2 — {vm2 < 244/^2)- For the system con- 
sidered (the values of the elastic moduli are given below) 
the second case is realized. At u < Vm2 Eq- (H) for a = 2 
gives 



y^±{\±iip), 



(7) 



where 



99 = y^(V^ + 62)72. 



(8) 



For the medium 1 the solution of Eq. (0) which satisfy 
the boundary conditions corresponds to the cases, for 
which Eq. (0) has two real and two imaginary roots or 
four imaginary roots 



j/1,2 = ±K = ± Y V -Di - 61 
2/3,4 = ±ii = ±iJ ^/D'^ + 61, 



(9) 



where Di = h\ — ci. (The k is a real value at v > C44/P1, 
and an imaginary value at v,ni < v < C44/ pi] u,„i is the 
root of equation Di[v) — 0) 

There are two types of the solutions of the wave equa- 
tions, which satisfy the boundary conditions (we will refer 
on them as I and II) 



ul{z) = Cni{z) 
ul{z) = ^C^/f (z). 



(10) 



where C^ is the normalization factor, p=I,II. In Eq.dlQ) 
fl.{z) is the odd function and fl{z) is the even function, 
while f^{z) is the even function and /"(z) is the odd 
function. 

The dispersion equation for v is obtained from the com- 
mon linear problem on the factors A'^^.. For the mode I 
the dispersion equation reads as 



i^itanh^ga taii^qa + R2t'An\\Kqa + i^stan^ga 

For the mode II it is modified to 
i?iCothK(7a cot^ga — i?2C0thK(7a + R:iCot^qa - 



i?4 = 0. 

(11) 



Ra 



= 0, 

(12) 



where 



Ri = {kitIk. + ^m^){ipm\ — Xm^) 

i?2 = m^m^{X^ + Lp^) - (p^iml + ml) 

-Kmf,{^m^ - (pm^) 

R3 = -m^m^iX'^ + ip^) + (puiml, + ml) 

+^m^{ipm^ - nm^) 

i?4 = (^m^ - ^mi^){ipm\ + Xm^) 

k(ci2 + C44) 

"^« = 2 \ 2 

CllK^ — C44 + PlW 
i{ci2 + C44) 



m^ 



-Cii^2 _ ^^^ ^_ p^y2 



m\ 



R = 



XciiiR- 1) 

C12 + C44 

_ ipC44{R+l) 

C12 -I- C44 
' cii{p2v'^ - c'li) 
C44{P2V'^ - C44)' 



Numerical solutions of Eqs. ( pl] , p^ ) for the parameters 
cii = 12.3, C12 — 5.7, C44 = 6.0 (all in 10^^ dyn7cm^) ver- 
sus the ratio between the thickness of the central layer 
(d = 2a) and the wavelength are shown in Fig. ||. One can 
see from the dependences presented, there are two gap- 
less modes (one of them is of the type I, while the other 
is of the type II) in the system. When the wavelength be- 
comes shorter, additional modes with higher frequencies 
emerge. The interaction with the gapless modes is only 
considered below. The quantities ff{z) can be written 
through elementary functions of z, but the expressions 
contain a complicate dependence on q and v. There- 
fore we will not present here the analytical expressions 



for ff{z). At an example, the plots of Ux{z), Uz{z) at 
d/l = 1 are shown in Fig. || . 




FIG. 1. Dependences of the velocities of the elastic modes 
on the ratio between the thickness of the central layer d and 
the wavelength I. Solid lines indicates the type I modes; 
dashed lines, the type II modes. 
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FIG. 2. Dependences of the displacement vector compo- 
nents (in relative units) on z (in units of the wavelength) at 
djl = 1; a, the type I mode; b, the type II mode. 



II. PIEZOELECTRIC COUPLING CONSTANTS 

Let us calculate the renormalization the velocity of the 
nonuniform elastic mode, caused by the interaction with 
the double-layer electron system with the coordinates of 
the electron layers zi(2) = ±a. We use the approach simi- 
lar to Ref . . Let us write the Hamiltonian for the elastic 
waves in terms of the phonon creation and annihilation 
operators (5+, h) 



H^ = ^tjg(6+6q 



^). 



(13) 



The Hamiltonian of the electron-phonon interaction is 
chosen in the form 



^int = -7= ^ / d'^r gqm'^rm.'^rme" 



qr 



(6q + 6lq), (14) 



where \I'+(\E') are the electron creation (annihilation) op- 
erators, TO, the number of the electron layer, 5*, the area 
of the layer. 

To obtain the matrix elements gqm we write the inter- 
action of the elastic wave with the electrons in the form 



H 



^^ , d^r et^rm^in^rm, 



(15) 



where iprm is the scalar potential of the electric field, gen- 
erated by the elastic wave in the layer m. The value of 
ip is determined by the solution of the Poisson equation 



Av3 = -{■iTr/e)l3ijkdiUjk, 



(16) 



where e is the dielectric constant, AjTc is the piezoelec- 
tric tensor, Ujk is the strain tensor. Under the choice of 
the X, y and z axes along the [100], [010] and [001] di- 
rections correspondingly the /3 tensor has nonzero (and 
all equal to the same value^) components ioi i y^ j ^ k. 
Under substitution of Eq. (|o|) into Eq. ( [l6| ) we obtain the 
following equation for the Fourier-component of the elec- 
tric potential 



where 



(9,2-g2)<^,(z) = z(47r/3A)CVW 



/(z) = q'fP{z) - 2qd.mz). 



(17) 



(18) 



We assume that the e and (3 parameters are same for 
the whole system. Then the boundary conditions on 
if reduce to the requirement of continuity of ifniz) and 
dz^q{z) at the interfaces. The solution of Eq. ( p7| ) with 
the boundary conditions has the form 



ipqiz) = i{A7rp/e)CPxpi^) 



(19) 



where xi{z) is an even function and xii(^) is an odd 
function. We do not present the analytical expressions 



for the Xp{z) functions here. As an example, the depen- 
dences fq{z) at d/l = 1 are shown in Fig. ^ (we use the 
values of C, calculated below) 



Dig, u) = (/ - i)(0) (g, cu)V (q))"' &-^^ {q, to), (25) 



where 




FIG. 3. Dependence of the Fourier component of the elec- 
tric potential (in relative units) on z (in units of the wave- 
length) at d/l = 1; 1, the type I mode; 2, the type II mode. 

To find the normalization factor we rewrite the com- 
ponents of the displacement vector in the form 



w,(r,z,t)-^CP/f(^)e"i--(6q-ffo^ 



(20) 



and substitute the expression (EG) into the common ex- 
pression for the elastic energy. Comparing the result with 
Eq. dl3), we obtain the following expression for the fac- 
tors CP 



where 



jP 



CP = {SvP)-^/'^, 



Aqj dzp(z)(|/P(z)p 
/o 



l/fWH- 



(21) 



(22) 



We note, that the q dependence of /p ( and C^ corre- 
spondingly) is reduced to the dependence on the param- 
eter qa. 

Rewriting Eq. ( p9| ) in terms of b operators and substi- 
tuting the result into Eq.(|l5|), we found the expression 
for the matrix elements in Eq. (HJ) 



gqm = ^{A^TPe/e^/vIP)xp(zr, 



(23) 



The velocity renormalization Atj and the absorption 
coefficient F for the elastic wave coupled to the two- 
dimensional electrons (the interaction is described by 
Eq.(p^) are given by the equation 






■i<l,^q)9m'q^ 



(24) 



where D^^^i is the electron density-density response 
function. The values of D in the random phase approxi- 
mation can be found from the equation 



VI 



■(9) 



2^2 



■(^W+(l-'5w)e-''^) - (26) 



the Fourier components of the Coulomb interaction, Z)^"^ 
is the density-density response function for the system 
without the Coulomb interaction. The quantities D'"' 
can be expressed through the longitudinal conductivity 
of the electrons 



i?^°^(9,-) 

mm ■ ' 



_ ^ rnm 

9 XX 

cue 



(q,^), 



(27) 



where a^^ = cr^^ are the diagonal with respect to the elec- 
tron layers, a^^ = cr^^, nondiagonal with respect to the 
layers components of the conductivity. Here and below 
we consider the case of two equivalent electron layers. We 
should note, that the nondiagonal components of conduc- 
tivity are equal to zero for the electron gas in the random 
phase approximation, while they may be nonzero for the 
composite fermion gas due to the interlayer statistical 
interaction (this case is considered in the next section) 
Substituting Eqs. (||, Hf^) into Eq. iM), we obtain 



ia+^{q,vq)/a 



Av F 
V q l + ialx[q,vq)/aj,j 



-a_- 



-icrxJq,vq)/aM 
1 + icrxx{q,vq)/a]^,j ' 



(28) 



where cr^^ = aH ± aH, a^j = ue/27r(l ± exp(-gd)) 



a± 



4nl3^ \Xp{a)±Xp{-aW 
ev-^lP 1 ± exp{—qd) 



(29) 



The a± functions play the role of the piezoelectric 
coupling constants, which are introduced for the con- 
sideration of the interaction of the SAW with the two- 



dimensional electrons. One can see from Eq. (29), that 
a+ coefficient is nonzero for the type I mode only, while 
the a_ coefficient is nonzero for the type II mode only. 

The dependences of a^ for the I mode and a_ for the 
II mode versus the parameter d/l are shown in Fig. 0. 
The parameters P = 4.5 • 10^ dyn^^^/cm, e=12.5 are 
used. The dependence of «_ on d/l for the transverse 
mode considered in Ref. ||] is also shown in Fig. ^. One 
can see from the dependences presented, the interaction 
of the electrons with the elastic modes, elliptically po- 
larized in the sagittal plane is much stronger then the 
interaction with the transverse mode. Note, that the 
case considered in has the advantage, that there is 
only one transverse nonuniform mode (if the thickness of 
the central layer is not very large comparing to the wave- 
length) , and its frequency is lower than the frequencies of 
the bulk modes. (For the case considered here the bulk 



transverse mode polarized along the [110] direction has 
the frequency, which is lower then the frequencies of the 
I and II modes. 



o 




FIG. 4. Dependence of the piezoelectric coupling constant 
on the parameter d/l\ 1, a+ for the type I mode; 2, q_ for 
the type II mode; 3, a- for the transverse mode [5]. 



III. ELASTIC MODE VELOCITY SHIFT UNDER 

THE PHASE TRANSITION IN THE 

FRACTIONAL QUANTUM HALL SYSTEM 

Let us apply the results obtained in the previous sec- 
tion to the study of a possibility of the observation of 
phase transitions in the double-layer fractional quantum 
Hall systems. To describe the quantum Hall system we 
use the composite fermion approach [the Chern-Simons 
fermionic model, developed in Ref. g for the double- 
layer system). Within such an approach the fractional 
quantum Hall system is modelled as a gas of the com- 
posite fermi-particles which carry the auxiliary statis- 
tical charge and the flux of the statistical gauge field. 
For the double-layer model two types of the statistical 
charges corresponding to two layer and two types of the 
gauge fields are introduced. In general case the composite 
quasiparticles carry the fluxes of the both types, namely 
the even number ^ of the flux quanta of the statistical 
field corresponding to their statistical charges and the in- 
teger number s of the flux quanta of the field, which cor- 
responds to the statistical charges in the opposite layer. 
This model at s y^ corresponds the states described 
by the Halperin wave function [p|. We refer the phase 
with s y^ as the phase with the interlayer statistical 
interaction. 

In the average field approxinration the partial screen- 
ing of the external magnetic field B emerges due to the 
influence of the statistical fields 



-Bn 



B{l-y{ij + s)), 



(30) 



where v is the filling factor calculated for one layer. The 
fractional quantum Hall effect corresponds to the filling 



factors, for which the value of BcH corresponds to the 
integer number TV of the filled Landau levels: 



N 



N{'ii) + s)±r 



(31) 



where the upper sign describes the case of i?off > 0, and 
the lower sign - the case of i?off < 0. One can see from Eq. 
(pi|), that certain fixed filling factors may correspond to 
different sets of the ■0 and s parameters (which describe 
different phases). 

To analyze the incompressible states (which corre- 
spond to the filling factors (|3l|)) it is convenient to ex- 
press the quantities a^x through the polarization tensor 
components H: 



^+(-) = _An+(- 






(32) 



where H+(-) = H" iH^^ H", H^^ are the diagonal and 
nondiagonal with respect to the layers components of the 



polarization tensor. The values of Hqq 
T+(-) _ e2g2 So 



+(-) 



H 



00 



27rwe A+(-) ' 



are found to be 

(33) 



where 



•^o — So 



nib 



^.^ (So(S.+TV)-S?), 
A+(-) = (l-(V'±s)Ei)2 

\2v^ /v^ Ar\ '^* ~ "^fe 



-(V'±s)^I]o(S2+iV) 



m*N 



(34) 

(35) 
(36) 



n 



xe 



N-l oo 

n=0 m=N ^ ' ^' 

'- m — n I 



S, = (Sgn(i?eff))^ 

a;™-"-i(m-n) 



m — n)^ 
dL™-"(x),, 



-2a;- 



dx 



(37) 



In Eqs.( |33| - |37l ) ujc = 2Tino / m* N is the effective cyclotron 
frequency, x — (qAoff)^/2, where Aoff = (A^/27rno)"^/^ is 
the effective magnetic length, L™^"(x) is the generalized 
Laguerre polynomial, m* is the effective mass of the com- 
posite fermions, mh is the band mass of the electrons, no 
is the average electron density. We used the modified 
random phase approximation [nO| for the calculation of 

n+(-) 



Substituting Eqs.(g2|,g3D into Eq.(|28D we find 



Au 

V 



E+So 



Eg Sq 



Eq Sq 



En So 



(38) 



r+(-) - 



where Eg ^~' = (e'^(7/ewc)(l ± exp{—qd)) (the absorp- 
tion coefficient F is equal to zero for the incompressible 

states). 



Let us consider the filling factor v = 1/5. If there is 
no interlayer interaction, this filling factor corresponds 
to the parameters ip = A, 5 = 0, N ^ 1. When the inter- 
layer distance becomes smaller a transition to the phase 
^ = 2, s = 2, A^=l may take place. One can see from 
Eqs. ( P4[]37| ), the first term in Eq. (|3q ) is not changed 
under such a transition. Therefore, the value of Av for 
the type I mode remains unchanged. On the contrary the 
jump of the phase velocity takes place for the mode of the 
type II (in this case the dependence on A_ which is the 
function oi ip — s survives in Eq. (BSh) The jump can be 
observed under the interlayer distance variation, but it 
can be hardly realized in experiments. But the effect can 
be observed indirectly if one measures the dependence of 
Av on the wave vector. If one use the elastic mode ve- 
locity at i/ = 1 as the bare value (at j/ = 1 the bare shift 
is determined by the same Eq. (|3^) with the parameters 
m* = rrib, 4' = ^ = 0, N = 1), the dependences of the 
relative shift on the wave vector distinguish qualitatively 
for the cases s = and s ^ 0. 

Fig. H illustrates this behavior. The dependences Av/v 
on the inverse wavelength (where Av is the difference be- 
tween the velocities at i^ = 1/5 and v = 1) are shown. 
We use the parameters uq = 10^^ cm^^, d — 500A, 
nib = 0.07r7ie, m* ~ Amb (m^ is the electron mass) for 
the calculations. 




FIG. 5. Dependence of the velocity shift on the inverse 
wavelength for the type II mode at the filling factor u — 1/5; 
1, the phase ^p = 4, s — 0; 2, the phase tp = 2, a = 2. 

One can see from Eq.(|3l|), that the phase transition at 
certain fixed filling factors may be accompanied by the 
change of the sign of i?off- For instance, dX v = 2/7 the 
phase without the interlayer statistical interaction corre- 
sponds to the parameters -0 = 4, s = 0, iV = 2 {B^s < 0), 
while the phase with the interlayer statistical interaction 
may correspond to the parameters ip — 2, s — 1, N — 2 
{Bes > 0). The transition between these phases results 
in a jump of the velocity for the modes of the both types. 
But the qualitative behavior of Av{q) remains practically 
unchanged. 

The mode of the type I, for which the maximum of 



the coupling is shifted to the long wave region, can be 
used for the indirect observation of the dependence of 
the effective magnetic length on the filling factor. The 
dependences Av{d/l) (the value of Aw is calculated rela- 
tive the velocity for the mode with the same I sA v — \) 
are shown in Fig. ^ at i^ = 1/3, 2/5, 3/7 {ip = 2 and 
N = 1,2,3 correspondingly). Here specify the case of 
s — Q. One can see from Fig. |[ the value of Aw oscillates 
as the function of l/l. The period of the oscillations is 
reduced when v approaches to the 1/2 value. It reflects 
that the effective magnetic length is increased for this 
sequence of the filling factors. 




FIG. 6. Dependence of the velocity shift on the inverse 
wavelength for the type I mode; 1, v = 1/3; 2, v = 2/5; 3, 
V = 3/7. 

Thus, the interaction of the double-layer electron sys- 
tem which is realized at the interfaces of the wide quan- 
tum well in the AlGaAs — GaAs — AlGaAs heterostruc- 
ture with nonuniform elastic modes, localized in the cen- 
tral layer of the heterostructure and elliptically polarized 
in the sagittal plane, is studied theoretically. The de- 
pendence of the piezoelectric coupling constant on the 
ratio between the thickness of the quantum well and the 
wavelength is found. It is shown, that the coupling con- 
stant increases under the decreasing of the wavelength 
and it reaches the maximal value at the wavelength which 
is of order of thickness of the GaAs layer (the concrete 
value of this wavelength is determined by the type of the 
nonuniform mode). The effect considered can be used 
for an experimental study of the dynamical properties 
of two-dimensional electron layers at large wave vectors, 
for which the interaction of the electrons with the sur- 
face acoustic wave is suppressed exponentially due to the 
finite value of the distance between the electron layer 
and the surface of the sample. The renormalization of 
the phase velocity of the nonuniform elastic modes cou- 
pled to the double-layer fractional quantum Hall system 
is calculated. It is shown, that the transition of the Hall 
system into the state described by the Halperin wave 
function results in that the dependence of the velocity 
shift on the wave vector is modified quantitatively for 



certain nonuniform elastic modes. 
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